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ABSTRACT 

We study the excitation of fine structure levels of C i, C ii and O i by ultraviolet (UV) 
photons around strong UV sources which are also ionizing sources of the cosmologi- 
cal reionization at redshift of ~ 10. The evolutions of ionized regions around a point 
source are calculated by solving rate equations for non-equilibrium chemistry. Signals 
of UV photons through the fine structure lines are considered to be stronger at lo- 
cations of more abundant chemical species of Ci, C ii and O I. Such environments 
would be realized where strong fluxes of non-ionizing UV line photons available for 
the pumping up of fine structure levels exist, and simultaneously ionizing UV photons 
are effectively shielded by dense H i regions. Signals from H i regions of moderately 
large densities induced by redshifted UV photons emitted at the point sources are 
found to be dominantly large over those of others. We discuss the detectability of the 
signals, and show that signals from idealized environments will be possibly detected 
by radio observations with next-generation arrays to come after the Atacama Large 
Millimeter /submillimeter Array (ALMA). 

Key words: atomic processes - hydrodynamics - ISM: Hii regions - cosmology: 
observations - dark ages, reionization, first stars - radio lines: general. 



1 INTRODUCTION 

At the epoch of Big Bang Nucleosynthesis of cosmic tem- 
perature r > 0.1 MeV, light elements of ^H,^'''He and '^Li 
are produced in abundances larger than 10"^" times that 
of hydrogen. As the temperature decreases, light elements 
recombine with electrons, and trace amounts of chemical 
molecules form during the Dark Ages between the recombi- 
nation and the reionization of the universe (Vonlanthen et 
al. 2009). The universe is then considered to experience the 
reionization at temperature of ~ 30 K, i.e., cosmic redshift 
of 2: ~ 10 (Larson et al. 2010). Sources of the reionization 
are perhaps ultraviolet (UV) photons from quasi-stellar ob- 
jects (QSOs), young galaxies and hypothetical low metallic- 
ity objects in the early universe, i.e.. Population III stars, or 
some exotic energy injection processes triggered, for exam- 
ple, by the particle decay. If stars including Population III 
have contributed to the reionization, the metal enrichment 
in the universe would occur simultaneously. The reioniza- 
tion history as well as the metal enrichment of the universe 
is, however, not yet determined precisely. Many suggestions 
have been made for investigations of the cosmic reioniza- 
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tion using atomic physics with future observations in radio 
frequency. 

The first method is observations of intergalactic 
medium (IGM) through the redshifted 21 cm line of neu- 
tral hydrogen (Madau, Meiksin & Rees 1997). The neutral 
gas regions could be quickly preheated to temperature larger 
than that of cosmic background radiation (GBR) in the pres- 
ence of a strong flux of UV photons through the photoion- 
ization heating (e.g., Chen & Miralda-Escude 2004). In such 
environments, the spin temperature determined by the pop- 
ulation fractions of the upper and lower states of hyperfine 
levels is forced to deviate from the CBR temperature by 
level-mixing via Lya photon scattering, i.e., Wouthuysen- 
Field effect (Wouthuysen 1952; Field 1958). There would 
be brief chances of observing them in absorption, and long 
chances in emission against the CBR. Future radio observa- 
tions at meter wavelengths of IGM would produce plenty of 
information on the epoch, picture and sources of the reion- 
ization. 

Although the 21 cm signal of the reionization of the uni- 
verse is expected to be relatively strong, the strong contam- 
ination from Galactic and extragalactic foregrounds makes 
it difficult to observe the absolute signal over the whole sky 
except for fluctuations around it. Measurements of angular 
fluctuations in the Hi signal are very difficult, while fiuc- 
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tuations in frequency would be easily separated and be de- 
tectable with future large arrays (Gricdin & Shaver 2004) . As 
for angular fluctuation, Di Mattco, Ciardi & Miniati (2004) 
have shown that angular fluctuations in 21 cm emission at 
^ > 1' can be detected if efficient subtraction of foreground 
sources is performed successfully. The 21 cm signal not by 
Wouthuysen-Fiold effect but by an effect of collision in ear- 
lier universe of higher densities has also been suggested to 
provide a valuable informations in the epoch when the spin 
temperature traced the gas kinetic temperature (Loeb & 
Zaldarriaga 2004). 

Long-lived inctastablc 26' state of hydrogen has been 
studied, and signals of the reionization through 3 cm fine 
structure line (2si/2 2^3/2) have been found to be unde- 
tectable with existing radio telescopes (Dijkstra et al. 2008). 
The reason is that a Ly/3 photon which is the pumping source 
of the ground state can excite only one hydrogen atom into 
the 2s state, while a Lya photon can excite many atoms 
into the excited spin state. The signals through hyperfine 
transitions of Dl, ^Hell (Dcguchi & Watson 1985) also ex- 
ist. Sigurdson & Furlanetto (2006) have suggested a chal- 
lenging possibility of determination of D /H abundance ratio 
by looking at hyperfine structure line at high redshift of 
both the H and D. Signals from IGMs and dense objects 
at high rcdsliifts through the "^He hyperfine transition has 
been studied. Detailed physics and suggested observational 
plans were presented (McQuinn & Switzer 2009), and the 
anticorrelation of H I 21 cm maps and ^He II 3 cm maps was 
suggested (Bagla & Loeb 2009). 

The infrared fine structure lines of Cl, C 11 and O I (Bah- 
call & Wolf 1968) play an important role in cooling of 
neutral regions (HoUenbach & Tielens 1999). Astronom- 
ical observations of those lines provide rich information 
about physical conditions in observed objects (Kaufman 
et al. 1999). The signatures of existence of heavy chemi- 
cal species like C and O through resonant scatterings with 
GBR imprinted on the GBR spectrum have been inves- 
tigated (Basu, Hernandez-Monteagudo &: Sunyaev 2004; 
Hernandez-Monteagudo, Verde & Jimenez 2006). With the 
Planck ^ spacecraft, the signatures can be searched for in the 
temperature anisotropy of GBR, and constraints on abun- 
dances of heavy species will be derived from the observations 
in reverse. 

Idea of UV pumping through fine structure lines of met- 
als during reionization as a way to distort the GMB dur- 
ing a reionization epoch has been suggested (Hernandez- 
Monteagudo et al. 2007) . They focused on the O I fine struc- 
ture, and provided estimations for the order of magnitude 
of the distortion for different levels of UV fluxes. The an- 
gular power spectrum of the distortion introduced by the 
Oi hosting regions at high redshift has subsequently been 
studied (Hernandez-Monteagudo et al. 2008). 

Observations of bright objects in radio frequency have 
been performed in search of emission lines of fine structures 
of heavy chemical species. The C 11 157/im line emissions of 
high redshift QSOs [J1148-h5251 at z = 6.42 (Maiofino et 
al. 2005), BR 1202-0725 at z = 4.7 (lono et al. 2006)] and 
a galaxy [BRI 0952-0115 at z = 4.43 (Maiolino et al. 2009)] 
have been detected. Walter et al. (2009) observed the high 
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redshift quasars (J1148-f5251 at z = 6.42) and derived an 
upper limit on a N 11 205/im line luminosity. The targets 
of these observations were, in fact, star- forming interstel- 
lar medium which emits through atomic fine structure lines. 
There are dense systems detected which absorb photon spec- 
tra of QSOs at fine structure transitions of C I, G 11 and their 
physical conditions have been constrained (Silva &: Viegas 
2002; Quast, Baade & Reimers 2002; Reimers et al. 2005). 

In this paper, we focus on cosmological signals of only 
UV sources without star forming activities through fine 
structure lines and their detectability. The structure of this 
paper is as follows. In Section 2. the pumping mechanism 
of fine structure levels and the dependence of the pumping 
efficiency as a redshift are shown. In Section 3, the model 
for calculation of reionization around a point source is in- 
troduced. In Section 4, the result of the calculation is pre- 
sented. Intensities of signals of UV radiation from the source 
through fine structure lines are estimated, and the depen- 
dence of signals on the redshift is shown. We discuss the 
detectability of signals, and claim that signals may be de- 
tected by future radio measurements in 0(100 — 1000 GHz) 
range. In Section 5, we summarize our conclusions. 



2 ULTRAVIOLET PHOTON PUMPING TO C I, 
C II AND O I METASTABLE STATES 

2.1 Populations of Fine Structure Levels 

2.1.1 Formulation 

For a pair of two states, i.e., a stable and an unstable states, 
the spin temperature Ts is defined as 



where rii {i — 1 for the upper state, i = for the lower) is the 
number density of the state i, gi = 2Si -|- 1 is the spin degrees 
of freedom of state i with spin Si, ft is the Planck's constant, 
1^10 is the frequency corresponding to the transition from the 
state 1 to 0, fe is the Boltzmann's constant. At the second 
equality, the temperature T* = huio/k was defined. The spin 
temperature describes the ratio of occupation degree of the 
upper state to that of the lower. 

For a group of three states, i.e., a ground state, a first 
excited state unstable to the decay into the ground state, 
and a second excited state unstable to the decay only into 
the first excited state, the spin temperature Ts2 is defined 
as 




where the subscripts i = 1 and 2 are attached to physical 
values of the first and second excited states, respectively, 

the temperature T,2 = hu-zi/k is defined with the frequency 
for the transition from the state 2 to 1, i.e., ^'21. Note that 
the spin temperatures for second excited states are defined 
using the number density of the i = 1 and 2 states, not i = 
and 2 states. 

Fine structure levels of Cl, C 11 and O I can be excited 
by continuum radiation at the frequency corresponding to 
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the energy levels of excited states, UV radiation at the fre- 
quency corresponding to the energy levels of intermediate 
bound states, and by collisions. For two discrete fine struc- 
ture energy levels denoted by and 1, the rate equation for 
their population is written as 



duQ 



drii 
~dt 



c^Jb gi 



2hv'^ go 



no 



2hvi 



■ niPi(^ + noCoi 



■ niCio, (3) 



where Aio is the Einstein ^-coefficient of the 1 — >■ transi- 
tion, 1/ is the frequency of the transition and do and Coi 
are the coUisional rates of the 1 — > and — i> 1 transitions, 
respectively. Jb = / 4>{v)IvdvdQ./ {^tv) is the mean intensity 
of radiation, where 0(f) is the line profile function and Ii, is 
the specific intensity at the transition frequency. The rates 
of excitation and deexcitation by UV photons, i.e., Pqi^ and 
Pqi^ are given by 



E 



AjoJjogj \ Aji 



50 



2hu 



AjiJjigj 
.91 



A 



3() 



(4) 



(5) 



where j indicates excited levels above the fine structure mul- 
tiplets and i is either or 1. There are several or more ex- 
cited levels of C I, C ii and O I where transitions from the 
ground state fine structure levels by electric dipole interac- 
tion are allowed. Although many lines could contribute to 
populating excited states of fine structures, we treat only 
transitions corresponding to the lowest energies for respec- 
tive species for the moment. All of the permitted lines are 
considered in Section 4.2.2. 

Figure 1 shows the transitions from the ground state 
fine structures to the lowest energy levels which is connected 
to the fine structures through strong electric dipole transi- 
tions. Numbers attached to lines are relative strengths of 
spontaneous emission rates A. 

The spin temperature in the steady state is derived 
from equation (3). We assumed the Planck distribution 
of the mean intensity Jb with the temperature of cos- 
mic background radiation (CBR), i.e., Tcbr = 2.725(1 + 
z) K (Mather ct al. 1999) with the rcdshift z. In addition, 
the thermal energy distribution of gas is assumed. The col- 
lisional transition rates are then given by 



Cij — Titargefyij (-^gas), 



(6) 



where ^-target IS the number density of target species and 7^ 
is the coUisional (de)excitation coefficient as a function of 

the gas temperature, i.e., Tgas. 

For the C 11 two level states, the spin temperature is 



Ts 



In 



J/LC 



+ yce 



where Tuv is given by the relation 
^ = ^exp^^" 



(tuv) 



(8) 



(Field 1958; Deguchi & Watson 1985), The temperature 
Teas satisfies 



Pio 90 

^01 yi \ -tci 



and 



pUV 

[1 - exp(-T./TcBR)] , 



A 



10 



(9) 

(10) 
(11) 



yc = ^[l- exp(-T,/TcBR)l 
were defined. 

For throe level states, the equations describing the bal- 
ance between production and destruction of the first and 

second excited states arc given by 



noPoi + W2-P21 = ni (Pio + P12) : 
noPo2 -I- niPi2 = 712 {P20 + P21) ■ 



(12) 
(13) 



where Pij includes transitions i — > j by all of the CBR con- 
tinuum, UV radiation and particle collision. The sum of the 

above two equations is equal to the balance equation for the 
ground state. These equations are solved for number frac- 
tions of the energy states, that is 

P01P20 + P21P02 -I- P01P21 



m 
no 

n2 
no 



P10P20 + P12P20 + P10P21 ' 

P10P02 + P12P02 + P01P12 
P10P20 + P12P20 + P10P21 ■ 



(14) 



(15) 



Spin temperatures, i.e., Ts and Ts2 of the first and second 
excited states are calculated using equations (1) and (2). 

Table 1 shows data on the fine structure and UV tran- 
sitions taken from Ralchenko et al. (2008) as well as data 
on the coUisional deexcitaion coefficients taken from Hol- 
Icnbach & McKee (1989). As for Cl the difl'crcnce between 
energy levels of ^P2 and ^Pi in the excited states higher 
than that of the ground state by ~ 6 x 10'' cm~' is 40.51 
cm-i (Ralchenko et al. 2008). 



2.1.2 Effects of UV pumpmg and Collisions 

We compare effects of UV pumping and collisions on popu- 
lations of fine structure levels. The efi^ects are not so large as 
that of CBR if the physical condition which is considered is 
not dense and not very near to extremely luminous objects. 
First, we define fij as the number ratio 



/ij(j/L, Tuv, J/c, Tgas, Tcbr) = Ui/rij 



(16) 



Using equation (7) we can describe the deviation of fu) for 
a two level system under a UV source and/or a coUisional 
environment from that for no such effects by the first order 
Taylor expansions in yi, and yc, i.e., 

A/io(yL, Tuv, yc, Tgas, Tcbr) _ /io(j/l, Tuv, t/c, Tgas, Tcbr) 



/io(0,^ Tuv, 0,^ Tgas, Tcbr 
■ C^) r / T* . T, 



/io(0,^ Tuv, 0,^ Tgas, Tcbr) 



2/L 



exp 



(- — 

V Tuv 



+yc 



exp - 



Tcbr 

n T 



- 1 



Tgas Tcbr 
(A/io//io)l + (A/io//io)c. 



- 1 



The first and second terms in the left hand side are con- 
tributions of UV pumping and collisions, respectively. The 
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Figure 1. Diagrams for ultraviolet permitted transitions which mix the ground state fine structures (solid lines) and do not (dashed 
lines). The attached numbers are relative strengths of spontaneous emission rates A. 



Table 1. Fine structure and UV transitions 

Fine structure transition 
Species State u (GHz) A (s"^) 



7H (cm3 s-i)° 



State 



UV transition 
u (cm"'^) 



A (s-i) 



Ci 


3P1 


-^Po 


492 


7.88 X 10-8 


1.6 X 10" 


-lOjiO.W 


3Pl 


-3Po 


6.035 X 10^ 


1.13 X 10* 






-3Po 


1301 




9.2 X 10" 


-11 j.0.26 


3Pl 


-3Pl 




8.64 X 10'' 






-3P1 


809 


2.65 X 10-'^ 


2.9 X 10- 


-lOj'0.26 


3Pl 


-3P2 




1.44 X 10* 
















■^P2 


-■^Pl 




8.58 X 10'' 


















_3p2 




2.52 X IQS 


Cii 


'P3/2 


-^Pl/2 


1901 


2.30 X 10-6 


8.0 X 10- 


-lojio.or 


'D3/2 


-^Pl/2 


7.493 X 10* 


2.37 X 108 












^D3/2 


-'P3/2 




4.74 X 10'^ 


Oi 




-3P2 


4745 


8.91 X 10-5 


9.2 X 10" 


-112^0.67 


^Sl- 


-3P2 


7.679 X 10* 


3.41 X 10^ 




3Po 


-3P2 


6805 




4.3 X 10- 


-lljiO.SO 


3Sl- 


-3Pl 




2.03 X 108 




3Po 


-3Pl 


2060 


1.75 X 10-5 


1.1 X 10- 


-107^0.44 
^2 


3Sl- 


-3Po 




6.76 X 10^ 



^ Deexcitation rate coefficients for collisions with Hi taken from HoUenbach & McKee (1989). T2 
rgas/(102 K). 



explicit expression for /lo value is 

/lo(yL, TuV, yc, Tgas, Tcbr) 
_ £1 r„-T./TcBR+yj^(e-T./T„v 



go 



+yc (' 



(18) 



Figure 2 shows contours of (A/io//io)l (upper panel) 
on the parameter plane of UV color temperature Tuv and 



scattering rate Pi/ for the UV transition P-- 



3/2 



D 



3/2, 



and also contours of (A/io//io)c on the parameter plane 

of the gas temperature Tgas and the number density of H I 
(lower panel). For this figure, the redshift was assumed 
to be 2: = 10, and the target in coUisional (de)excitation 
was assumed to be H i. Solid and dashed lines correspond 
to positive and negative values for indicated amplitudes, re- 
spectively. Because of Tuv S> T, the effect of UV pumping 
is not sensitive to Tuv [cf equation (17)]. (A/io//io)c is 
positive in the parameter region of Tgas > Tcbr while it is 
negative in the parameter region of Tgas < Tcbr [cf. equa- 
tion (17)]. 

Figure 3 shows contours for C I fine structure lines of 
(A/io//i())l (sohd lines) and (A/2i//2i)l (dotted) (upper 
panel) on the parameter plane of Tuv and Pi, for the UV 
transition ^Po (ground state) — Pi, and also contours of 
(A/io//io)c [solid (for positive values) and dashed lines 
(for negative values)] and (A/2i//2i)c [dotted (positive) 
and dot-dashed lines (negative)] on the parameter plane of 



Tgas and nui (lower panel). The redshift was assumed to 
be z — 10, and the target in coUisional (de)excitation was 
assumed to be H I. 

Figure 4 shows contours for O I fine structure lines of 
(A/io//io)l (solid lines) and (A/2i//2i)l (dotted) (upper 
panel) on the parameter plane of Tuv and for the UV 
transition '''P2 — Si, and also contours of (A/io//io)c 
[solid (for positive values) and dashed lines (for negative 
values)] and (A/2i//2i)c [dotted (positive) and dot-dashed 
lines (negative)] on the parameter plane of Tgas and tihi 
(lower panel). The redshift was assumed to be = 10, and 
the target in coUisional (de)excitation was assumed to be 
Hi. 

From Figs. 2, 3 and 4, it is found that under environ- 
ments of intense UV fluxes and relatively low densities and 
low gas temperatures considered in this study, the effect of 
UV pumping is typically larger than that of collisions. We, 
therefore, neglect the effect of collisions on the population 
of fine structures in what follows. 



2.2 Intensities of Redshifted Fine Structure 
Radiation 

We estimate a signal from a region including abundant C i 
or C II, or O I of angular size on the sky that is larger than a 

beam width, and of radial velocity range that is larger than 
the bandwidth. For a fine structure line of frequency fo, a 



© 20XX RAS, MNRAS 000, 1-23 



Signals of the cosmological reionization in the radio sky through C and O fine structure lines 5 



10 



10 



of 10 



10 



10 



-14 

1000 





- 


2x10"= 


- 


r 


2x10"* 






2x10"^ 






2x10"" 




2x10 



10' 

Tuv (K) 



10" 



10 



0.1 



s 0.01 



10 



10 



1 ' 1 

: / 
/ 1 




; 

: / 

: / 1 




; 
/ 


l^^^^^^^^^^^^^xKT^^ 


^ ^ 1 




1 

/ 


\ \ 2x10"* 


2x10"'' 






\ \ 2x10"^ 


1 


\ 2x10"*^ 

i . ^ . 1 1 



10 



100 
. (K) 



1000 



10* 



Figure 2. Contours of deviations in number ratio ni/ng from 
the values for only CBR contribution. The upper panel shows 
the contribution of UV pumping on the parameter plane of UV 
color temperature Tuv ^^'^ scattering rate for the UV tran- 
sition ^P3/2 ^3/2- The lower panel shows the contribution 
of collisions on the parameter plane of the gas temperature Tgas 
and the number density of H I njjj . Solid and dashed lines corre- 
spond to positive and negative values for indicated amplitudes, 
respectively. 



intergalactic optical depth at frequency vo/{l + z) along the 
line of sight is (Hernandez-Monteagudo et al. 2007) 



Aij\%nA„ (z) (gifj/gj - fi) 



(19) 



where {i,j) = (1, 0) or (2, 1), Xij is the wavelength for transi- 
tion i j, nA„ (z) is the number density of chemical species 
An (C I, C II or Oi) at redshift z, fi = Ui/ (X^ji-i) is the 
population fraction of energy state i, Ho = Tl.Olj'^ km 
Mpc~^ (Dunkley et al. 2009) is the Hubble constant, Q^n is 
the energy density parameters of matter (e.g., Madau et al. 
1997; Padmanabhan 1993). The optical depth in the homo- 
geneous universe is typically much less than unity. In the 
optically thin case, the intensity of the signal AJ^ with re- 
spect to that of the CBR radiation B^{Tcb'r.) is 



Ah 



exp(T,/rcBR) - 1 



exp(r./r5) - 1 



T, 



(20) 
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Figure 3. Same as in Fig 2 but for n\/no [solid (for positive 
values) and dashed lines (for negative values)] and n2/n\ [dotted 
(positive) and dot-dashed lines (negative)] of C I. The upper panel 
shows the contribution of UV pumping on the parameter plane 
of Tuv ^iid Pf for the UV transition ^Pq (ground state) — Pi. 
The lower panel shows the contribution of collisions on the plane 
of Tgas and tihj • 



which is identical to equation (9) in Hernandez-Monteagudo 
et al. (2007) under the assumption that the background radi- 
ation field is that of thermal CBR. In this study, we investi- 
gate signals of the UV continuum radiation left on the radio 
background via excitations of fine structure lines. The UV 
source is assumed to have a power law luminosity spectrum, 
i.e., Lv <x . When the UV photon scattering occurs 

close to the emission redshift, the total scattering rate Pi, 
for a point source is given by 



>^gu L„/{hv) 
OTT gi Anr-' 



'Aui\ ("L. 



'a, 47 



(r-Mpc/0.1) 



(cf. Madau et al. 1997) where the subscripts I and u 
mean the low energy state of fine structure and the up- 
per energy state to which the transition from the state I 
is allowed, respectively, r is the proper distance between 
the emission and scattering. In the second equality, \a ~ 
1.215 X 10"^ cm is the wavelength for Lya transition of H I, 
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Figure 4. Same as in Fig 3 but for O I. The upper panel sliows 
the contribution of UV pumping on the parameter plane of Tuv 
and Pn for the UV transition ''P2 Si - The lower panel shows 
the contribution of collisions on the piano of Tgas and njjj - 



Ac = 6.265 X 10* is the spontaneous emission rate for 
Lya, (i/L^)^ is the product of the luminosity of the point 
source and the frequency of photon at Lya frequency in units 
of 10*'' ergs s~^, and rjupc is the distance r in Mpc. 

Figure 5 shows signals through fine structure lines rela- 
tive to corresponding cosmic background radiation as a func- 
tion of the redshift under the following assumptions: The 
universe is homogeneous, and the number density of hydro- 
gen is nn = 1-9 x 10~^ cm~^(l -I- z)'^ . The abundances of C 
and O in the universe at redshift 2 are t/c = 1-4 x 10~* 
and yo = 3.2 x 10~*, respectively, which correspond to 
local ISM values (Maiolino et al. 2005)^. The elements 
(C and O) are assumed to be completely in the chemical 
species (Ci, Cii and Oi), i.e., ncj — nnyc, = nnyc 
and noi = nnyo, respectively. A QSO of luminosity of 
L„ — 4.1 x 10^^(i//i/a)~"^ ergs s^^ Hz~^ [corresponding 
to (i/L^)^^^ — 1] with as = 1/2 (solid lines), as = 3/2 



^ Although the abundances of C and O are thought to be smaller 
in the earlier epoch of the universe, their chemical evolution as 
a function of the redshift is rather uncertain. We, then, adopt 
here the present values for the abundances which would be upper 
limits on the old day abundances. 




-10 



Figure 5. Signals through fine structure lines in the homoge- 
neous universe relative to corresponding cosmic background ra- 
diation as a function of the redshift. A QSO of luminosity of 
= 4.1 x 10''^'^{u/i'a)'°'S ergs s^^ Hz^i with as = 1/2 (solid 
lines), 03 = 3/2 (dashed) and as = 5/2 (dotted) is assumed to 
shine at a distance of 1 Mpc. The respective lines are derived un- 
der the assumption that the elements (C and O) are completely 
in the chemical species (C I, C II and O l), and therefore maximum 
values. 



(dashed) and as = 5/2 (dotted) is assumed to shine at a 
distance of 1 Mpc (rMpc = 1). A decrease in flux of UV 
photons associated with continuous scattering with C I, C 11 
and O I during their propagation is not considered. This ef- 
fect is found very important in Section 4 below. In drawing 
this figure we confirmed that the spin temperature of O I 
228 K line calculated with equations given in Sec 2.1 is ex- 
actly equal to that described by an approximate equation 
(7) of Hernandez- Monteagudo et al. (2007). 

At high redshift of z > 10, all signals are small at the 
distance of rjvipc = 1- As the redshift decreases, spin tem- 
peratures become smaller with the CBR temperature de- 
creasing by redshift. Optical depths decrease since the num- 
ber density of particles decrease with the cosmic expansion 
[see equation (19)]. Since the number density of the CBR is 
smaller and its energy is lower at low redshift, the UV pump- 
ing effect gets relatively dominant over the CBR effect below 
some critical redshifts for respective fine structure lines. Spin 
temperatures then decouple from the CBR temperature at 
the critical redshifts in the present setting of fixed radius. 
Since the CBR temperature at low redshift is less than the 
line temperatures, i.e., (T* > Tcbr), the decoupling of the 
spin temperature (Ts > Tcbr) leads to the large difference 
of line intensity from that of CBR. The factor in the angle 
bracket in the right hand side of equation (20) is, therefore, 
large. The signal to CBR ratios are enhanced at low redshift 
for these reasons. 



2.3 UV Color Temperature 

We mention an important point on the UV color tempera- 
ture, for example, of the ground state of C I. 

In a time scale for UV photons to scatter C I particles 
by UV pumping, i.e., the mean free time, the color tempera- 
ture of the UV photon around the line center relaxes to the 
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kinetic temperature of the gas (i.e., Ci here) (Field 1959; 
Rybicki & dcU' Antonio 1994). The mean free time of a line 
center photon is t''""^ = l/(nciCr^r) with aj? the UV line ab- 
sorption cross section at the line center. The line absorption 
cross section is 

a, = ^ijLA^icj>{u). (22) 
877 gi 

Neglecting the microturbulence and the very small effect of 
natural broadening, the line center cross section is 



,2 

^ul 9u < 



91 ■^/ttIS.i/u 



where 



c V m 



(23) 



(24) 



is the Doppler width with m the mass of an atom (Rybicki 
& Lightman 1979). 

Using this equation and nci = n^ShUc the mean free 
time is found to be 

^■"• = «„^£,(i±i)-'c(jj|^)-'. (25) 

where Tgas,K is the gas temperature in units of K, 5h is the 
local density excess in units of the universal average value 
of baryonic matter. The color temperature around the line 
frequency, therefore, could approach to the gas temperature 
depending on the temperature of the gas, redshift, density 
excess, and the abundance of the scattered species. Since the 
Doppler width Ai^d = 2.2 x '^Q^s^^T^i^^ is much smaller 
than the frequencies corresponding to the hyperfine transi- 
tions oiv = 0(10^^ — lO^^s"^), line photons of different fre- 
quencies are completely distinguished at scattering processes 
with C and O. In other words, a photon with frequency for 
one line transition can not induce scatterings corresponding 
to other transitions. Furthermore, strong lines which enable 
changes in photon energy do not exist around the frequencies 
of the UV pumping lines. The color temperature is, there- 
fore, given by an original color temperature of the source 
Tuv = hi\ine/[k(i + as)\ before relaxing to the spin temper- 
ature by scattering with C or O. This situation is completely 
different from the case of the Lya pumping in which Lya 
photon can mix the hyperfine levels of the H I ground state 
many times (Madau et al. 1997). 



3 NUMERICAL CALCULATION OF ID 
HYDRODYNAMICS 

We calculated a hydrodynamics coupled to a chemical reac- 
tion network assuming the spherical symmetry of space in 
order to show an example of geometric structure of ionized 
region. 

3.1 Hydrodynamics 

3.1.1 Fluid Equations 

The hydrodynamic equations to be solved in this study are 
as follows: 



|+V.(pv) = 0, 



|^ + (v.V)v = -^-V$, 
at p 



d{pE) 

at 



+ V-[{pE + p)v]=S, 



(27) 



(28) 



where t and r are the time and the position, respectively, 
p, V, p and $ are the density, velocity, pressure and grav- 
itational potential as a function of t and r, respectively, 
E = v'^ /2 + e is the fluid energy per unit mass with 
e = p/ [(7 — l)p] the internal energy and 7 the ratio of spe- 
cific heat, 5 = r — A is the energy gain with F and A the 
heating and cooling terms, respectively. The ratio of specific 
heat is fixed to be 7 = 5/3 in this study. 

We here introduce the comoving coordinates (Shandarin 
1980; Shapiro & Struck-Marcell 1985) using comoving vari- 
ables to describe deviations of the fiuid quantities from those 
of the expanding homogeneous universe. New variable de- 
noted by primes are 



dt' = a ^dt, 


(29) 


1 -1 
r = a r. 


(30) 


/ 3 
p = a p. 


(31) 


v' = a [v - H{t)v\ , 


(32) 


p = cC'p, 


(33) 


e' = o^e = p / [(7 - 1) p'] , 


(34) 


S' = a'^S, 


(35) 



where a{t) — 1/(1-1-2;) is the scale factor of the universe with 
z the redshift, and H{t) — a /a is the Hubble expansion rate. 
The peculiar gravitational potential is assumed to be negli- 
gible, and only the potential in a homogeneous isotropically 
expanding universe is considered here, i.e.. 



$(t,r) 



2%Gp{t)r^ 



(36) 



where G is the gravitational constant, and p{t) is the mean 
density at time t. 

We assume that the cosmic expansion is described 
by ACDM model. The cosmological parameters arc taken 
from a fit to the Wilkinson Microwave Anisotropy Probe 
(WMAP) 5 year data (Model ACDM-hSZ-Flens). (Dunkley 
et al. 2009)^^ The energy density parameters of baryon, 
matter and dark energy are Qh ~ 0.0441 ± 0.0030, Sim = 

0. 258 ± 0.030 and = 0.742 ± 0.0030 respectively. In this 
study we concentrate on the early epoch of the universe, 

1. e., z > 1. The model, therefore, is effectively the CDM 
dominated universe. 

Under the assumption of the CDM dominated universe 
and of the ratio of specific heat 7 = 5/3, the fluid equations 
are transformed into the following ones: 



i;+v'.(pv)=o. 



d{p'^') 
dt' 



+ V 



\pv 



+ P 



0, 



(37) 

(38) 



(26) 



^ http://lambda.gsfc.nasa.gov. 

Adopted parameters are almost the same as those of WMAP 7 
year data (Larson et al. 2010). 
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dt 



P' 



= S'. 



(39) 



where E' = v' /2 + e' was defined. These equations do not 
include the cosmic expansion explicitly . They are solved 
simultaneously with the time evolution of the scale factor 
a{t). 



3.1.2 Heating and Cooling Rates 

As a heating term, the photoionization is included. The heat- 
ing rate from kinetic energy of photoelectrons and protons 
produced by photoionization is given by 



r = 



^ Ui X 

=Hi,Hei,Heii 



/>oo 



dE^ {E^ - Eth,i) n^{E^)ar''{E^)c, 



(40) 



where is the energy of nonthermal photon, Eth,i is the 
threshold energy for photoionization of species i, n-y{E.y) 
is the energy spectrum of the photon number density as a 
function of time t and position r, af"^{E^,) is the photoion- 
ization cross section, and c is the light speed. The production 
of secondary electrons at the photoionization of H and He is 
followed, and additional energy inputs by these electrons are 
taken into account in the estimation of the heating rate. For 
the photoionization cross sections we used the FORTRAN 
subroutine phfit2 by D. A. Verner^. We adopted the fitted 
functions from ShuU & van Steenberg (1985) for heat genersr 
tion firactions associated with secondary electrons produced 
by the H I and He I photoionization. 

As cooling terms, radiative coolings due to coUisional re- 
combination, excitation, ionization and free-free transitions 
as well as the Compton cooling off the GBR are included. We 
adopt the cooling rates for Case B (optically thick case) re- 
combinations of Hll and Hem by Hummer (1994) and that 
of Hen by Hummer & Storey (1998). The cooling from the 
dielectronic recombination of He ii is also included with its 
rate taken from Shapiro & Kang (1987). The cooling rates 
for the collisional line excitation are calculated using the ex- 
citation rates given in Sobelman, Vainshtein & Yukov (1981) 
for excitations of Hi, Hel and Hen from the ground states 
of main quantum number n = 1 to excited states of n = 2,3 
and 4. 

Below 10* K, metastable transitions or fine structure 
transitions of dominant metal species contribute to the cool- 
ing. We, therefore, include important line transitions of C I, 
Cii and Oi (Hollenbach & McKee 1989). The cooling rate 
for fine structure transitions of species i is given by 



Afe = mj^ fi^ [eM-Eui/Tg^) - exp(-£;„,/rcBR)] 

gi 



Ku 

i=Hi,e 



(41) 



where 7^; is the deexcitation rate of transition -u — > Z at colli- 
sions with species j whose value was taken from Hollenbach 



'' Subroutines by D. A. Vcrncr are 
htt p : // www . pa. uky. edu/ " verner/fortran. html . 



available at 



& McKee (1989). For this equation, the detailed balance be- 
tween forward and backward reactions is assumed, and the 
spin states of species are assumed to be described by the 
CBR temperature. (Note that spontaneous emission rates, 
i.e., Ay,i are typically laxger than collisional rates, i.e., 7^; 
in low density environments considered in this paper.) The 
cooling rate for (semi) forbidden transitions of species i with 
transition temperatures of 0{W^ K) is given by 

Afor = V ^ [exp(-£„o/rgas) - exp(-Ko/TcBR)] 

^-^ go 



(42) 



3=Hi,e 



where it is assumed that almost all particles of species 
i are in its ground state, i.e., /o = 1, and the term 
exp(— _Euo/Tcbr) is negligible. These metal transitions are, 
however, found not to affect this calculation significantly. 
The cooling rate for the collisional ionization is given 

by 



Aio 



E 

=H,,Ho,,Hc„ 



ili{(riv) 



col 



(43) 



where Ii is the ionization energy of species i, and {atv)''"^ is 
the thermal average value of the cross section times velocity 
for collisional ionization. The {aiV^"^ values are calculated 
with a Verner's subroutine, i.e., cfit. The cooling rates for 
the free-free transitions and the Compton scattering of CBR 
by free electrons are taken from Shapiro & Kang (1987). 



3.1.3 Numerical Method 

We use the Monotone Upstream-centered Schemes for Con- 
servation Laws (MUSCL) with the Roe's method and the 
second order Runge-Kutta method in time integration. The 
computation is then second order accurate in space and 
time. The number of grid point is 400 and the spacing is 
Ar = 0.5 kpc at the redshift of « = 8.7. The computational 
region is thus ^ r ^ 0.2 Mpc. 

The calculation domain is assumed to be uniform in 
density initially. A point source is put on the origin and 
starts lighting at t = 0, whose luminosity is given by 

= 4.1 X 10^1 ergs s"^ Hz'^ (y^Ybv) ' 

where as is the spectral index fixed to be 1.5, and 10.2 cV 
is the energy corresponding to the Lya transition of H I 
atom. The amplitude and spectral index of the point source 
are the same as those in Madau et al. (1997). This point 
source emits total ionizing photon for hydrogen of = 



13.6 eV/h 



L^/(hv) dv = 1.6 X lO'^** S" 



3.2 Chemical Reaction Network 

We constructed a chemical reaction network code including 
22 chemical species and 55 chemical reactions. The species 
are atoms and all ions of H, He, C and O as well as the 
electron. 

The rate equation which was solved is given by 

dm ^ /%\ ^ /dyi_\ ^ /dyi_\ ^ / dyi_\ ^ 

dt \ dt ) ion \ dt J rec \ dt J col \ dt J cha ' 
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where yi = rn/riH is the relative number density of the 2/0 = 3.2x10 are used, which correspond to the local ISM 



ith species to that of Inrdrogcn nu = nui + nun, and the 
right hand side includes terms for the nonthermal photoion- 
izations (ion), thermal radiative recombinations (rec), and 
thermal coUisional ionizations (col) of all atoms and ions, 
and charge transfer reactions (cha) (e.g. Kingdon & Ferland 
1999). Each term is described by 



(dyi\ 
V dt /io 



yj=(i+e-) f dE^nj{E^)ariEj)c 



POO 

-Vi / 



dE^n^{E^)a'r(E^)c, (46) 



dy 



'j =(i+e-) 

dt 



= nnyiVeicrivY'"' 
-nnyjy^icTjvy 



(dyi\ I \c 



-niiyiye{crivY° , 



-nnyiyk{<yikjvT ^) 



(47) 



(48) 



(49) 



where j = (i + e~) is the species which is photodisintcgratcd 
into the species i and an electron, {aivY'"^ is the thermal 
average value of the cross section times velocity for recom- 
bination. The first and second terms in the right hand sides 
of the above equations correspond to the production and 
destruction terms, respectively. In the fourth equation aijk 
is the cross section of reaction i + j ^ k + I with any I, and 
indexes j and k should be summed. 

The Verner's subroutine rrfit is used to evaluate the 
radiative recombination rates. The collisional ionization by 
secondary electrons produced in the photoionization of H 
and He is taken into account. We adopted the fitted func- 
tions from ShuU & van Steenberg (1985) for the ratios of 
energies used for secondary ionizations to the total energy 
produced by the photoionization of Hi and He I. Charge 
transfer rate coefficients arc taken from Kingdon & Ferland 
(1996, and updates and addenda)*^: FORTRAN subroutines 
by Jim Kingdon are used to calculate reaction rates of H 
and all ions of He and C. Fitted coefficients by Rakovic et 
al. (2001, unpublished) for reactions of H and all ions of O 
are adopted. Charge transfer reactions with He of C"*" with 
q = 3 and 4, and C"*" with q = 2 — 4 are also included with 
their rate coefficients taken from the web page. 



3.3 Input Physics 

The initial abundances of H, He, C and O arc given as fol- 
lows. The helium to hydrogen number ratio is j/hc ~ 0.082 
as predicted in standard big-bang nucleosynthesis model 
with baryon-to-photon ratio determined from WMAP, i.e., 
7? = 6.3 X 10-1° (Dunkley et al. 2009). yc = 1.4 x 10"* and 



^ Updates and addenda are shown at http://www- 
cfadc.phy.ornI.gov/astro/ps/data/. 



values (Maiolino et al. 2005). All particles arc assumed to 
be in the form of neutral atoms. The initial number density 
is uniformly given by uh = 1.9 x 10-^5b(l + z)^ cia~^ cor- 
responding to the adopted f2b value. Two cases of 5h = 10^ 
and 10* are calculated and presented in this paper. We as- 
sume that there is no astrophysical heat source before the 
epoch corresponding to the beginning of the calculation. The 
initial gas temperature then should be given by that hav- 
ing decoupled from the CBR temperature at z ^ 200 and 
experienced the redshift as T oc {1 + z)^. We adopt such 
temperature calculated by Loeb & Zaldarriaga (2004), i.e., 
T = 2.3K [{l + z)/10f. 

The local ionization flux as a function of frequency z/ is 

F^{r) = L^. exp(-r^)— !--//(2;start - Zi), (50) 

where H{x) is the Heaviside's step function, Vi is the the 
frequency which redshifts to v after running a distance r. fi 
is, therefore, the initial frequency of photon emitted at the 
origin at redshift Zi, i.e., 

Vi 1 + Zi 



1 + z 



(51) 



8.7 is the initial redshift in the calculation, Tv(r) is 
the optical depth at time t between the origin and a given 
position r, and is properly given according to the position 
r. For example, when the radius r is larger than all the 
ionization fronts, it is 

r^(r) ~ nn [(r - rf') (yniihv) + ynei {r - rf"') a^eiihv) 



■rf'^")'^He"„(M], 



(52) 



where r\ is the radius for the ionizing front of species i, i.e., 
the boundary between the ionized and un-ionized regions. 
The first, second and third terms in the square bracket in 
the right hand side are for ionizations of H I in the region of 
r ^ r^' , He I in the region of r ^ r^*' , and He ii in the He ii 
region enclosed in r^'^" = = rf^^ , respectively. 

Dusts might have been produced by Type II supernovae 
even in the early universe (Dwek et al. 2007). Recent dis- 
coveries of hot-dust free quasars aX z ^ & may indicate pos- 
sible differences in dust properties in early epoch of the uni- 
verse (Jiang et al. 2006, 2010). The dust abundances in the 
early universe is, however, still very uncertain. The effects 
of dust are then not considered in this study. In the energy 
region of E~, ^ 7.5 cV (corresponding to the frequency of 
A ^ 0.17 yum) related to the UV pumping, an estimate of 
mean extinction Text of the diffuse ISM with solar metallic- 
ity indicates Text/iVn ^ 0(10"^^) cm^ (Zubko et al. 2004) 
where A'h is the hydrogen column density. If abundances 
of dust are similar to that in the solar neighborhood, the 
extinction length scale would be 



[Text/iVH] nu 



^ 0(0.1 Mpc) 



U03; 



(53) 



3.4 Grid Size 

In order to calculate evolutions of the H ll region precisely, a 
grid spacing should be small enough to correspond to an op- 
tically thin distance (Madau ct al. 1997). Madau et al. (1997) 
have set the ionizing photon flux to be zero for energies less 
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than 1.8 ryd, the He I ionization potential to prevent too 

rapid a growth of the Hi ionizing front. In the present set- 
ting, the distance corresponding to optical depth of Atv = 1 
is given by 

-3 



> 10- 



(54) 



where the ionization cross section of a^' ^ 6 x 10~^** cm^ 
[e.g., fig. 2.2 in Osterbrock (1989)] was used. The grid size in 
this calculation is ~ 0.5 kpc. This is optically thick for 5h = 
10"^ or 10*. Although optically thin grid spacings thus could 
not be used in this treatment of coupled hydrodynamics and 
chemistry, we do not focus on details of the evolution of 
ionizing front in this study. 

The approximated treatment taken here is to account 
for an optical depth in the photo-induced heating and reac- 
tion rates. We input the average number density of photon 
n-y{Ej) over the grid interval region into equations (40) and 
(46). The average number density over the grid interval be- 
tween r and r + Arigrid is given by 



jA.i,Hd F,{r)/{hu) exp [-r.(drO] dr' 
F^(r) 1 -exp(-Ar^) 



(55) 



h'^u At, 

where Atu is the optical depth between r and r + Arigria. 

3.5 Time Scales 

The time scale of light front move is 

The time scale of recombination in the fully-ionized state is 

Atrec = 



7 X 10^ Myr ri> «« (1 + 2j/He)-^ " S^' 



-13 T-i-0.845 



(57) 

cm^ s~^ (Madau et al. 



where as ~ 2.6 x 10 

1997) is the recombination coefficient to the excited states 
of hydrogen, and T4 is the gas temperature in units of 10* K. 

The ionization time scale at around the ionization front 
of Ti/ = is 

-1 

Atu 



nHi _ 
10 yr 



n.y{E-^)a'^^{E-y)c dE^ 



0.1 Mpc 



(58) 



3.6 Evolutions of Ionization Fronts 

The ionization time scale [equation (58)] is much smaller 
than the scale now treated in the hydrodynaniical calcu- 
lation. The photoionization, therefore, can be regarded to 
occur instantaneously. The total number of photon causing 
the ionization satisfies 

rt-At,{rl') 

VS'^ / fi{U)dU 
Jo 



An ( 



ri 



01 



ft frl'(t) 



A '2 0/ 

47rr neUjf^i^aB.i 



where S-y is the emission rate of ionizing photon assumed 

to be independent of time here, fi{t) is the number frac- 
tion of photon emitted at time t which is used for the 
ionization of species i, ae.i is the recombination coeffi- 
cient to excited states of species i, and two quantities, i.e.. 



rl/a and n°[iy 



=H for i =Hi and 



=He for i =Hel and Hell] , are measured in the co- 
moving frame. n^|.^j is the initial abundance of neutral Hi 
or He I. We note that for the case of He, a two-step ionizar 
tion of He l He 11 — > He lll need two ionizing photons per 
one He atom inside t he H ell ionization front. The quantity 
AU = [rj7(3c)(3i/()v/ar/2)^/^ + t^/^f - to is the time it 
takes for the light to travel the distance rl in the matter 
dominated universe with to the age of the universe at the 
initial time of the calculation. The left hand side of the equa- 
tion is the total number of photons used for ionization. The 
first term in the square bracket in the right hand side is the 
number of photons required to ionize the region of radius 
rl, while the second is number of times of recombination 
occurring inside rl. 

A time evolution of ionization front of species i, i.e., rl 
can be described by a solution to the equation of the total 
number of photon causing the ionization of i, i.e.. 



Jo 



where 



-At,(r;') 



fi^{t»)dt» 



47r {rl'Y 



Or 



(60) 



■5fk)Pk{i,t) 



^/fe/r(t)Pfc(i,t)(6i) 



is the effective fraction of photons emitted at t to ionize 
species i eventually, ft is the fraction of photons emitted 
in energy range k which is specified by Ei ^ E^ < E2 
(fe=l), E2 ^ E-f < E3 {k=2) and E3 ^ E-y {k=3), where 
El = 13.60 eV, E2 = 24.59 eV and E3 = 54.42 eV are 
the threshold energies of Hi, He I and Hell, respectively. 
5fk is the absorption factor by photoionization evaluated 
at ionization front of i. /I""" = 1 — Sfk/fk is therefore, the 
survival fraction of photons in energy range k. Pk{i,t) is 
the fraction of photons emitted at time t to ionize species i 
which is calculated by a prescription given in Appendix A. 

The ionization cross sections have peaks just above 
threshold energies and rapidly decrease as the energy in- 
creases [equations (Al), (A2) and (A3)]. We, therefore, ap- 
proximately identify recombinations of Hll, Hell and Hem 
as losses of photons in energy ranges of fc = 1, 2 and 3, re- 
spectively, as a result of subsequent ionizations by photons 
in the energy ranges. The amplitude of photon fiux in en- 
ergy range k at the ionization front is reduced by a factor 
fk^'^i't) in the present treatment. The recombination of each 
species then involves photons in only one energy range, and 
a relation holds: 



fnt) = 1- 



fkS-f 



rl'it) 



dr' ^4nr' 



2 0/ N 

nenj(^iyaB,i I 



dtt. 



d{t* - At*i) ' 



where i =Hi (for k = 1), Hei (for fe = 2) and Hell (for 
A: = 3). When all photons of E\ ^ E-, < E2 are used for 
recombination inside the Hi ionization front (/f""^ < is 
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then derived from naive estimate) , photons of E2 ^ Ej < E3 
are supplcmentarily used with /i"' fixed to be zero. The 
ratio dtt/d{tt — Afti) derives from the time integral range 
of the left hand side and the second term of the right hand 
side of equation (59). We evaluate ionization fronts at any 
t. They are determined only by ionizing photons arriving at 
the ionization fronts by the time t which should be emitted 
hy t- Ati. 

As for the input for fk"^{t) values, we assume that the 
temperature of ionization regions is Tgas = 3 X 10" K which 
corresponds to the pressure of P ^ 1 eV cm^ in Hii and 
Hem regions of 5b ~ 10^ (see Fig. 8). The number density 
of electron is approximately given by 



= (1 + 2yH. 

(1 + J/He)nH 
J/HeriH 




(0 ^ r < r 
Hen < 



Hcj 

I 



in"' S r) 



r < r 
r < r 



Using the fixed temperature and the electron number den- 
sity given above, /^"'(t) is calculated. 

Equation (60) is rewritten in the form of 



3& 



3 



(z) Jo 



t-Ati(rj') 



ff{U) dU = 0. 



(63) 



The integral is numerically estimated discretely, i.e., 

t-Ati(ri') 



ff{U)dU ^ ^lff{U-AUi)5{U-AUi)U64) 



where ff^{t — AU) is the effective fraction for photons emit- 
ted at time t — Ati, which arrive at the ionization front of 
i at time t. rl{t) values are obtained by solving equation 
(63). fi^''{t) values are simultaneously obtained. The solu- 
tion rl{t) is used as the evolution of the ionization front in 
this calculation. 

Figure 6 shows the radii of the ioinzation fronts for H I 
(solid lines), Hei (dashed) and Ifeii (dotted) as well as the 
radius of the Ught front night- Curves for three cases, i.e., 
5b = 1, 10^ and 10^ are drawn. The curves for 5b = 1 al- 
most overlap that for the light front. This means that there 
is only a narrow region where the hydrogen exists as a neu- 
tral atom inside the light radius of ~ 0(0. 1) Mpc. In this 
case signals of UV radiation through C I, C 11 and O I would 
be very weak since their abundances are very small after a 
very short time interval between the UV photon arrival and 
the ionization. We then calculate the evolution of the chem- 
ical abundances in the cases of 5b = 10^ and 10*, in which 
relatively large regions of neutral hydrogen exist inside the 
light radius since it takes long time to ionize dense regions. 
The density excess values of 5b — 10"^ and 10* seem some- 
what high. The cases, however, would quite roughly simulate 
situations where high density clouds exist in the line of sight 
towards the UV emitting object. 

Figure 7 shows the calculated survival fractions of pho- 
tons in three energy ranges after traveling from a QSO to 
ionization fronts, i.e., fk"^{t) for cases of 5b = 10'^ and 10*. 
The three energy ranges are as follows: (1) Ei ^ E^ < E2, 
(2) E2^E-y < E3 and (3) E3 ^ E^. The survival fractions 
are used in this calculation to take account of reduced UV 
fluxes at ionization fronts. 
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Figure 6. Radius of the light front and those of the ionization 
fronts for H I (solid lines) , He I (dashed) and He II (dotted) . Lines 
correspond to the cases of density excess values of 5b = 1> 10^ 
and 10*. 
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Figure 7. Survival fractions of photons in energy ranges k after 
traveling from a QSO to ionization fronts, i.e., /^"'^(i) in the cases 
of density excess values of 5b = 10^ a^nd 10*. 



4 RESULTS 

4.1 Evolution of an Ionized Region 

^.1.1 Case o/5b = 10^ 

Figure 8 shows the pressure (upper panel) and the velocity 
(lower) as a function of the radius at time t = 0.12 Myr 

(1), 0.24 Myr (2), 0.35 Myr (3), 0.47 Myr (4) and 0.59 Myr 
(5). A bump in the pressure curve appears at r ^ 0.015 
Mpc although it is difficult to read it clearly. This is caused 
by the discrete treatment of space. The point corresponds to 
the innermost grid where a light front abundant in energetic 
UV photons arrives faster than ionization fronts. Environ- 
ments rich in energetic photons make it possible to heat up 
gas particles efficiently. Inside this bump, gas particles re- 
ceive UV photons both of relatively high and low energies si- 
multaneously leading to a relatively inefficient heating. The 
bump causes a discontinuity in velocity. Although these dis- 
continuities are unphysical, they do not affect the result of 
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Figure 8. Pressure (upper panel) and velocity (lower) as a func- 
tion of the radius r at time t=0.12 Myr (1), 0.24 Myr (2), 0.35 
Myr (3), 0.47 Myr (4) and 0.59 Myr (5) for the case of (5b = 10^. 
Solid (dotted) lines correspond to calculations with grid number 
400 (800). 

structures of ionized region at larger scales. We then do not 
search for a solution to the bump here. 

Dotted lines in Fig. 8 show the same physical quan- 
tities for another calculation with a larger grid number of 
800 which was performed for a shorter time. Differences in 
pressure, or temperature, of the two calculations with dif- 
ferent spatial resolutions are seen to be small. We then do 
not expect so large errors in line signals derived from the 
calculated temperatures, as described below. Nevertheless, 
more complete estimations with high resolution calculations 
are desired. 

Figure 9 shows chemical abundances of hydrogen, he- 
lium and electron (top panel), carbon (middle) and oxygen 
(bottom) as a function of the radius r a,t t — 0.35 Myr (3). 
Thin solid lines show the same physical quantities calculated 
with 800 grid points. 

4.1.2 Case o/ (5b = 10" 

Figure 10 shows the pressure (upper panel) and the velocity 
(lower) as a function of the radius at time t — 0.12 Myr (1), 
0.24 Myr (2), 0.35 Myr (3), 0.47 Myr (4) and 0.59 Myr (5). A 
step in pressure exist at r ~ 0.04 Mpc (the ionization front 



T 




r (Mpc) 

Figure 9. Chemical abundances of hydrogen, helium and electron 
(top panel), carbon (middle) and oxygen (bottom) as a function 
of the radius r at t = 0.35 Myr (3) for the case of 5b = 10'^ • 
Thick (thin) lines correspond to calculations with grid number 
400 (800). 



of Hell). Sharp decreases in pressure cause some structures 
in velocity which are seen in Fig. 10. Dotted lines show the 
same physical quantities for another calculation with grid 
number of 800. 

Figure 11 shows chemical abundances of hydrogen, he- 
lium and electron (top panel), carbon (middle) and oxygen 
(bottom) as a function of the radius r at t = 0.35 Myr (3). 
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Figure 10. Pressure (upper panel) and velocity (lower) as a func- 
tion of the radius r at time t=0.12 Myr (1), 0.24 Myr (2), 0.35 
Myr (3), 0.47 Myr (4) and 0.59 Myr (5) for the case of <5b = lO". 
Solid (dotted) lines correspond to calculations with grid number 
400 (800). 



Thin solid lines show the same physical quantities calculated 
with 800 grid points. 



4.1., 3 Effect of UV background 

A UV background can cause a heating of gas. In addi- 
tion, even a weak UV background in the epoch of the 
hydrogen reionization had quickly ionized neutral car- 
bon with low ionization potential. For example, con- 
sider a region of density excess (5b = 10" at 1 -f z = 



10 exposed to a soft UV background with flux J — 
10"^^J2i(-B-^/13.6 eV)-^ '^ ergs s"^ cm'^ Hz"^ str"^ in 
E-y ^ 13.6 eV. Assuming the gas temperature of 300 K and 
that the abundance of electron is equal to that of C 11, the 
equilibrium abundance ratio of C 11 to C I is approximately 
ncu/^Ci = 3.8 X 10'* J21. Even a J21 value as low as lO""*, 
which is possibly achieved very soon after the onset of Pop- 
ulation III star formation (see e.g. Haiman, Abel & Rees 
2000) is enough to ionize most of the carbon. In order to 
check effects of UV background the following extra calcula- 
tions are performed. 

We assume that there are UV sources which ap- 
pear 0.1H{z = 8.7)"^^ before the onset of the point 
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Figure 11. Chemical abundances of hydrogen, helium and elec- 
tron (top panel) , carbon (middle) and oxygen (bottom) as a func- 
tion of the radius r at 4 = 0.35 Myr (3) for the case of = 10*. 
Thick (thin) lines correspond to calculations with grid number 
400 (800). 



source. H{z) * is the Hubble time at redshift z, and 
the flat ACDM universe has been assumed, i.e., H{z) = 
H(y^ r2m(l + z)'^ + S7a. The UV background is, therefore, 
switched on at z = 9.8 corresponding to the cosmic age 
of 0.50 Gyr. 
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The UV spectrum is roughly given by 



10-^^21 (E-,/13.6 eV)-^-^ ergs/s/Hz/str 

for E^f S 13.6 eV 
10-21 J2iex(£7/13.6 oV)-^-'^ crgs/s/Hz/str 
X exp[-1022(<_. + O.OSajj- 

for Ej > 13.6 eV, 



(65) 



where cr]°^^2 is the photoionization cross section of i in units 

of cm^. This form is adopted from Haiman ct al. (2000) al- 
though the modulation factor is neglected, and the spectral 
index is different. The exponential factor corresponds to the 
absorption above ionization energy threshold of H I by the 
neutral IGM under the assumption of a hydrogen column 
density of A^h = lO'^'^ cm"'^. A parameter ex is the ratio of 
X-ray to UV flux: ex = for pure stellar sources and ex ~ 1 
for typical quasar spectra. The parameter ex would thus 
parametrize fractions of contributions of stars and quasars. 
Two cases of ex = (soft) and 1 (hard spectrum) with a 
fixed value of J21 = 1 are calculated assuming the homoge- 
neous UV backgrounds. 

In both cases of e = and 1, C I is completely ionized 
after the UV exposures over 0.1 Hubble time at z — 8.7 
(88 Myr). In the case of e = 1, the ionization by hard UV 
photons heat up the homogeneous gas. The temperature at 
z = 8.7 is then 10^ K. We first calculated one-zone chemical 
abundances in an exposure to the UV background from z = 
9.8 to 8.7. Wc then calculated chemical abundances around 
the point source as a function of radius taking account of 
ionizations by the point source and the homogeneous UV 
backgrounds using results of one zone chemical evolution as 
input. 

Figure 12 shows the pressure (upper panel) and the ve- 
locity (lower) as a function of the radius at time t = 0.12 Myr 
(1), 0.24 Myr (2), 0.35 Myr (3), 0.47 Myr (4) and 0.59 Myr 
(5). Solid, dashed and dotted lines correspond to cases of 
a soft, hard and no UV backgrounds which is the same as 
in Fig. 8, respcctivch'. Calculations for two cases with UV 
backgrounds are performed with 200 grid points. There is 
no large difference between values in three cases except for 
pressure values in the case of hard UV background. 

Figure 13 shows chemical abundances of hydrogen, he- 
lium and electron (top panel), carbon (middle) and oxygen 
(bottom) as a function of the radius r at t = 0.35 Myr (3). 
Solid, dashed and dotted lines correspond to cases of a soft, 
hard and no UV backgrounds which is the same as in Fig. 
9, respectively. In the cases with UV backgrounds, C I has 
already been ionized while in the case of no UV background 
C I exists outside the light front. 



4.2 Signals through Fine Structure Transitions 

We divide signals by the UV photon pumping into two 
classes. One is a signal emitted in regions just outside of 
ionized regions. The other is that of UV photons generated 
as more energetic than that of the line center frequency and 
redshifted to the line center when they arrive at radii larger 
than ionization fronts. 




Figure 12. Pressure (upper panel) and velocity (lower) as a func- 
tion of the radius r at time t=0.12 Myr (1), 0.24 Myr (2), 0.35 
Myr (3), 0.47 Myr (4) and 0.59 Myr (5) for the case of <5b = 10^. 
Solid, dashed and dotted lines correspond to cases of a soft, hard 
and no UV backgrounds, respectively. 



4-2.1 Signals from Ionization Fronts 

The UV photons can pump up the energy levels of C I, C 11 

and O I when they arrive in regions of abundant given chem- 
ical species. Since the amount of UV photon is finite, flux 
densities of the pumping photons are reduced as they pump 
up. Firstly, we check if there are enough UV photons for 
mixing the fine structure levels so that the UV pumping 
can effectively occur inside the light front before arrivals of 
ionization fronts in the setting of this study. 

The ionization time scale at around the ionization front 
(i.e., an optically thin region) is given by equation (58), and 
the precise number regarding Hi ionization is Ation(Hi) = 
2.2 X 10* s [r/(0.1 Mpc)]2. Minimum cross sections for scat- 
tering off of UV photons by C i, C 11 and O i required for a 
scattering in the time scale of ionization, i.e., Ation(Hi) is 

1 



= 8.7 X 10-1^ cm2f^)"7i±i)' 

\io-*J V 9.7 J 



^(loO Vo.lMpc 



(66) 
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Figure 13. Chemical abundances of hydrogen, hehum and elec- 
tron (top panel), carbon (middle) and oxygen (bottom) as a func- 
tion of the radius r at 4 = 0.35 Myr (3) for the case of = 10^. 
Solid, dashed and dotted lines correspond to cases of a soft, hard 
and no UV backgrounds, respectively. 



The frequency ranges around the line center of frequency 
I'VY which satisfy a ^ 10"^'' cm^ in the case of Tgas = 
2.3 K[(l -I- z)/10]'^ (Loeb & Zaldarriaga 2004) is Av/uw = 
6.2 X 10~^ for the ground state C I atom, 6.8 x 10~^ for the 
ground state Cll ion, and 5.2 x 10~^ for the ground state 
O I atom. The emission rates S'g of UV hne photons which 




r (Mpc) 

Figure 14. Differences between spin temperatures and the CBR 
temperature (K) as a function of the radius r (solid lines). The 
value is negative for the O I 99K line at r ^ 0.07 Mpc due to 
the large population of second excited state of O I by an efficient 
pumping, and its absolute value is showrn (dashed line). 

interact with C and O before ionizing photons could ionize 
Hi are derived by 

Ss = ^Au. (67) 

Values in the present setting are obtained using equation 
(44): Ss = 3.8 X 10^^ s"^ for Cl, 3.0 x 10^^ s^^ for Cll, and 
2.2 X 10^^ s"^ for Ol. 

Figure 14 shows differences between spin temperatures 
and the CBR temperature in units of K as a function of 
the radius r. In drawing this figure, losses of UV line pho- 
tons by interaction with atoms and ions are not considered. 
As will be discussed below, UV line photons are quickly 
used to pump up the fine structure levels. This figure, then, 
shows the maximum temperature differences adequate for 
only regions right beyond the ionizing front where the loss 
of photons can be neglected. 

Suppose that ionizing UV photons are mainly used to 
ionize H I, and that photons of frequency around the UV 
pumping lines are mainly used to excite and deexcite C I, C ll 
and O I. We compare the emission rate of ionizing photon per 
hydrogen number density, i.e., fniS-,/nu with that of line 
photon per differences of number densities of excited states 
i, i.e., Ss/An^^ . Here An\^ = n^^ (Ts) - n^i^ (Tcbr) was 
defined, and number densities of excited state i, (T) 
are described by the steady state spin temperature (Ts) and 
that of CBR (Tcbr), respectively [see equations (1) and (2)]. 
At r = 0.1 Mpc, 5s/Ancj = 9.2/Hi5'T/nH and Ss/An^^ = 
lAfHj^Sj/nn for Cl, Ss/Atiq^^ = 77fiiiSj/nH for Cll, and 
Ss/Anl^ = 311/Hi5.,/nH and Ss/An^^ = 248/HiS.,/nH for 
O I. It is seen that UV photons which can scatter C and O 
exist in abundances larger than ionizing photons when nor- 
malized to target number densities. The scattering regions, 
therefore, extend outward faster than ionized regions. 

Inside ionization fronts of C I, C ll and O I, UV line pho- 
tons propagate outward without scattering. When they en- 
ter in region with abundant C I, C ll and O I, however, they 
scatter the species and change the fine structure population. 
The UV line photons effectively pump low energy level states 
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up to high energy levels until steady states realize. Excited 
states eventually decay into low states by magnetic dipolc 
photon emission. Repeating pumpings and decays lead to 
net reduction of number flux of energetic UV line photon 
which can pump up C and O from low energy states. 

The radius scale for this scattering loss is given by those 
of regions where the steady state spin temperatures can be 
hold by the flux from the point source. Balances between 
pumping up and spontaneous decay of excited states in a 
volume V = 47rr| Ar, with tb the boundary of the C and O 
ionization and Ar the width, are roughly described by 



4nrt 



( 



-Att) — ArAn 



•i^Ai,) ^ 0, 



(68) 



\hiAjY 

where the UV flux in the frequency range of Doppler width 
Aj/d is considered in the source term which enter into this 
region from the boundary re- The radius scale is thus 



Ar ' 



Az^ 



For the radius of r = 0.1 Mpc, this equation yields 



(69) 



gas,] 

",3\-lTnl/2 



Ar = 0.8 {Sh/lQ-'y'T^!^^ pc (for the Ci 24 K line), 



0.02 [5^/WT^T^L\ pc (Ci 39 K), 

0.02 {5^./lQ'r^Tlil^ pc (Cn 91 K), 

4 X 10-3 (<5b/10')-'rgY,'K pc (Oi 228 K), 

2 (5b/103)-^T^a^'K PC (Oi 99 K). 

The UV line photons around the line center quickly scatter 

C and O and decreases during propagations through narrow 
regions if the density is high as in this calculation. 

4-2.2 Signals from Outside of Ionization Fronts 

(i) Steady state abundances of excited states 
Secondly, signals from excitations by redshifted UV pho- 
tons are estimated. For simplicity, we assume that UV pho- 
tons redshifted beyond certain critical frequencies from the 
blue side instantaneously scatter C and O, and are used for 
pumping up of the species. The number flux of the UV line 
photon, i.e., Fi/jjY / {^^^y^ ■) should satisfy 



dr \hv\jY J 



iAjvH{t) 



(70) 



where AnA„ is the abundance of excited state i, Nea is 
the effective number of lines which contribute to pump- 
ing up C I, C II and O I of fine structure levels, and de- 
fined below. The first term in the right hand side is for 
spontaneous emission of excited state i, while the second 
is for production by redshifted UV photons. The produc- 
tion term is proportional to the number flux of the UV 
line photons and the rate of redshift per unit distance, 
i.e., Idv/dr] « \di'/{cdt)\ = vuy H{t)/c. The steady state 
for abundances of fine structure levels would be realized 
in a region between an ionization boundary and the light 
front. They are described by a balance in equation (70), i.e., 
d[F^„^/{hiAjv)]/dr = 0, as 



vvvH{t) 



(71) 



(ii) Efi^ective number of available UV lines 
UV photons described by a color temperature Tuv repeat 
scattering C and O species through the UV photon pump- 
ing, and eventually relax to the color temperature of the 
fine structure spin states which is approximately given by 
the CBR temperature if the UV photon pumping occurs fre- 
quently enough. In this process, chemical species in excited 
states which have been pumped up by UV line photons can 
emit line photons corresponding to the fine structure tran- 
sitions. The effective number iVeff is given as follows taking 
account of the relaxation. 

(a) Cii 

For the C ll, the net increase in the number of ions in 
the excited state is given by the net increase in the number 



of photons corresponding to the transition P; 



3/2 



3/2 



(see Table 1) after the relaxation process through the UV 
scattering. A^efi is, therefore, given by 

2 



iVeff = 



l-Fexp(-T,/rcBR) 



-1, 



(72) 



where the first term in right hand side is the relative 
number of line in the final relaxed state, and the second 
term is that in the initial state, that is equal to unity. In 
the final state, initial UV photon fiux for two lines, i.e.. 



3/2 



D3/2 and 



1/2 



lines by the equilibrium fraction of 



D3/2 is distributed to the two 



3/2 



:exp(T./TcBR), 



(73) 



1/2 



which realizes the balance between pumping rates of both 
directions, i.e., 



noJofo->u->i = fii^iCi->u->o, 



(74) 



where Fq = Fip^^^, Fi = F2p^^^, and ai^u^j = 
a{i ^ u) X ^iij/(X]fc ^uk) is the cross section of pro- 
cess i ^ u ^ j [cf. equation (22)]. 
(b) Oi 

Similarly, effective numbers for the first and second ex- 
cited states of O I axe given by 



1 + exp (-T,2/Tcbr) [1 + exp (-T,/Tcbr)] 



-1,(75) 



eft, 2 



3exp(-T*2/TcBR) 



1 + exp (-T.a/TcBR) [1 + exp (-T./Tcbr)] 



-1,(76) 



respectively. Here it was assumed that the UV photon 
pumping of the first excited state of O l is efficient as well 
as that of the ground state as a sufficient condition for 
the equilibrium UV line photon. Since the abundance of 
the excited Ol ('^Pi) state is rather small [goni / (gino) ~ 
exp(-228 K/Tcbr) = 2.3 x lO-"* ior 1 + z = 10], this 
assumption is satisfied in environments of high O I den- 
sities. If the density excess is 5h ~ 10'', the equilibrium 
UV flux can realize [equation (82)], and a contamination 
from a coUisional signal is still small if the gas tempera- 
ture is not so high (Section 2.1.2). If the density excess is 
smaller, the number of UV lines through which UV pho- 
tons effectively scatter becomes smaller. If the density of 
O I is as small as (5b ~ 1, only the UV pumping of the 
ground state is efficient. Signals of the Ol fine structure 
lines then originate from pumping of the ground state by 
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a part of UV line photons. In such a case of small O I num- 
ber density (Hernandez- Monteagudo et al. 2007), signals 
of both O I 228 K and 91 K would become emissions and 
their intensities smaller than the estimations given here, 
(c) Ci 

As for C I, two contributions through the upper levels 
^Pi and ''P2 exist. The effective numbers for the first and 
second excited states are then given by 

3 



1 + exp ( 



+ 



-r*2/TcBR) [1 + exp (-T./Tcbr] 
2 



i-^exp(-r.2/rcBR) 



A^o£F.2 = 



3 exp (— T,2/rcBR) 



1 + exp (-r.2/TcBR) [1 + exp ( 
2exp(-T.2/rcBR) 



-n/Tc 



+ 



(77) 



(78) 



i + exp(-r.2/rcBR) 

respectively. The first and second lines in right hand sides 
of the above two equations correspond to the contribu- 
tions of UV pumping through the ^Pi and ■^P2 levels, 
respectively. Since the number density of excited state C I 
(=^Pi) is large [goni/igmo) ~ exp(-24 K/Tcbr) = 0.41 
for 1 -I- 2 = 10], the equilibrium abundances realize rela- 
tively easily. 

(iii) Scattering rate 

A steady state scattering rate [cf. equation (21) for the 
rate derived neglecting the scattering loss] in the region out- 
side of ionization front is estimated as follows: The balance 
between a rate of scattering and that of depletion of UV 
line photons whose frequencies are redshifted over critical 
frequencies is described as 



nA„P^ = iVcff- 



(79) 



hvuY c 

We obtain the steady state P^, value from this equation, i.e. 



P. 



-iVc 



47rr2 



= 3.3 X 10" 



A^cff 



1 + z 



10 



10- 



(rMpc/0.1)2 \ 3.7 X 10-1^ s- 



(80) 



See Section 2.1.2 for efi^ects of UV photons on spin temper- 
atures. 

(iv) Signals through the lowest-energy UV lines 
Figure 15 shows the differential antenna temperatures, 
i.e., the brightness temperature Tb minus CBR temperature 
in units of mK as a function of the radius. For this figure 
we assume that the abundances of C I, C 11 and O I are their 
maximum values, i.e., yci = 1.4 x lO"**, ycu = 1-4 x 10"*, 
and yoj = 3.2 x lO"'', respectively. The differential antenna 
temperatures are derived with the following equation which 
is satisfied when Tb - Tcbr < Tcbr: 



7b — 2cE 



Ah 



exp(r./TcBR) - 1 



(81) 



Tcbr B„{Tcbr) (T* /Tcbr) exp(r* /Tcbr 

Since the differential antenna temperature is rather small, 
signals of these magnitudes would not be detected with ex- 
isting radio telescope or the planned Atacama Large Mil- 
limeter/submillimeter Array (ALMA). 
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Figure 15. Differential antenna temperatures (mK) as a function 
of the radius. Solid curves are for positive values, while dashed 
curves are for absolute values of negative ones. The respective 
curves are derived under the assumption that the elements (C 
and O) are completely in the chemical species (Cl, Cll and Ol), 
and therefore indicate maximum values for respective lines. 



(v) Signals through all available UV lines 

So far only the effects of lines shown in Fig. 1 are analyzed. 
There are, however, many other lines of energies less than 
the ionization threshold of H I which can mix fine structure 
levels. We then include effects of those lines. Lines avail- 
able for pumping up the C and O species in the fine struc- 
ture levels should satisfy a requirement: the time scale of 
UV pumping of C or O by photons of line center frequen- 
cies Atpump is shorter than that of redshift At^ad- If this is 
not satisfied, only a part (~ At^cd / Atpump) of the UV line 
photons can scatter C or O before their frequencies are re- 
moved from the line center and they become inert. Using 
Afpump =^/TTAi/-D/{'nAN'^) (Rybicki & Lightman 1979) and 
Atrod = Au-D /[ui!vH{t)], and defining A by ct ~ X^/{8tt)A 
[cf. equation (22)], a relation is obtained, i.e. 



At 

pump 
Atrf,d 



3 



10- 



A 



(82) 



where I'uv.is is the line center frequency in units of 10^^ s~^. 
The lines contributing to the pumping up of fine structure 
levels should have Einstein yl-coefficient enough large in or- 
der to satisfy Atpump/Afred < 1- This condition tends to 
be met in environments of high ShYA^j values. If the two 
timescales are nearly identical, i.e., Atpump ^ At^ed, the 
depletion at the UV pumping of UV line photons balances 
with the production via cosmological redshift. In such a sit- 
uation, the steady state UV line flux and the resulting line 
scattering rate [equation (80)] are relatively large. 

Assuming that all UV lines of A > lO'^ s~^ contribute to 
the pumping up, for example, signals through the C and O 
species are calculated using atomic data from Ralchenko et 
al. (2008) . Figure 16 shows differential antenna temperatures 
when contributions from all UV lines of ^4 > 10^ s~^ are 
included. The signals for C I and O I are enhanced relative 
to curves in Fig. 15 since many additional strong lines exist. 
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r (Mpc) 

Figure 16. Same as in Fig 15 but including the effects of all UV 
lines of spontaneous emission rates larger than 10® s^^. 

On the other hand, the signal for C ll is not enhanced much 
since only two upper states, one of which is the ^^3/2 shown 
in Fig. 1, are available. 

4.3 Flux Densities of the Point Source 

Total flux densities, i.e., S — J du, from IGM via the 
fine structure line emission or absorption are estimated with 
the calculated results of the time evolution of the chemical 
structure. We suppose that a point source starts lighting at 
time t — at redshift z — 8.7. The comoving distance to the 
source, rs, is defined by 

r S± S_ r (83) 
Jts ^0 ^n^a + (1 - n^)ai 

where as and ts are the scale factor and the corresponding 
time of the universe when the source exists, and ta is the 
age of the present universe. 

We introduce two dimensional Cartesian coordinates. 
The source and the observer are assumed to locate at the 
origin and x = 00. The flux density in the redshifted fine 
structure frequency which is measured on the earth is given 
by 

(1- + Z) rg J-x„ax 

where ro 0.2 M pc is the ra dius of t he calculation domain, 
a;max =y/r^ — y^, and r =y/x^ + is the radius from the 
source. We assume that the UV lines with emission rates 
of ^ ^ 10^ s~^ contribute to the UV pumping of C or O 
when the abundance is enough large, i.e., Ya^ ^ 10""*. The 
value of An\^{r,t + x/c) is thus given by equation (71) if 
Vajv (^) t+x/c) is larger than 10~*. The value is set to be zero 
otherwise. Note that the effect of the time delay is included 
in An^jy {r, t + x/c) since ionization fronts propagate at sub- 
light speeds. 

Figure 17 shows the total flux densities emitted through 
the C II and O I fine structure lines in the calculated domain 
of r ^ 0.2 Mpc in units of njy GHz. The calculated result 
for the case of 5^, = 10^ is used for this estimation. As the 
light propagates outwards, the volume of the region experi- 
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Figure 17. Total flux densities (nJy GHz) emitted through the 
Cji and Oi fine structure lines and the Hi 21 cm line in the calcu- 
lated region for the case of 5i, = 10^ . The values which we would 
measure at 2 = are shown as a function of the time from the 
start of the hydrodynamical calculation. Solid lines correspond 
to emissions, while a dashed line is for absorption. Thick lines 
correspond to the standard case of no UV background, while thin 
lines correspond to two cases of a soft and hard UV backgrounds. 

encing the UV pumping of C and O increases. The emissions 
(solid lines) and the absorption (dashed line) triggered by 
redshifted UV photons are then enhanced. For this figure we 
assume that the chemical abundances as a function of the 
radius are fixed after the time of the end of calculation, i.e., 
tend at values of tend- This is because we do not have results 
of t > tend which are necessary to calculate the total flux 
densities due to the time delay effect in equation (84). No 
C I region exists inside the fight radius in our calculation. 
There are large amounts of ionizing photons of C I escaping 
from shielding since the threshold energy of Cl (11.3 eV) 
is lower than that of Hi (13.6 eV). Thick lines correspond 
to the standard case of no UV background, while thin lines 
correspond to two cases of a soft and hard UV backgrounds 
(Section 4.1.3). 

In Fig. 17, the total flux density for the Hi 21 cm line 
is also shown. The fiux density is given by 

^ nl ,4^''''^'' I "y^y I "^'^^ AnUr,t + x/c), (85) 

where 1/21 = 1.4 GHz is the frequency corresponding to the 
21 cm line transition, A21 — 2.9 x 10~^^ s~^ is the spon- 
taneous emission rate of the hyperfine structure transition, 
and ArzH is the difference of the number abundance of the 
excited level of Hi hyperfine structure (njj) in the presence 
of the UV radiation field from that without the UV field. 

Hi regions which are irradiated with UV photons are 
heated predominantly via photoionization. The heating rate 
is given by equation (40). The gas temperature in the Hi 
region with enough UV flux would, therefore, be higher while 
that in region with strongly shielded UV flux would be lower 
than the GBR temperature. Note that the Lya scattering 
provides a negligible heating (e.g., Chen & Miralda-Escude 
2004). 

The spin temperature of the H I hyperfine structure, i.e., 
Ts, is usually much larger than the transition temperature 



© 20XX RAS, MNRAS 000, 1-23 



Signals of the cosmological reionization in the radio sky through C and O fine structure lines 19 



T, = hi'21/k = 0.068 K. The abundance of the excited state 
is then given by 

1 3exp(-T,/rs) 3/ ir,\ 

= l + 3exp(^T,/rs) ""- " 4 - iTs) ^^^^ ^''^ 

We define Ts and Tbg as the spin temperature in the 
UV radiation fields and that without it. Altliough the ef- 
fect of the collisional excitation of the hyperfine structure 
line is not negligible in regions of high densities, we here ne- 
glect the effect in order to derive the maximum effect of the 
Lya pumping on the 21 cm emission. Under this assump- 
tion, when there is no UV radiation field, the spin tem- 
perature is purely described by the CBR temperature, i.e., 
Tbg ~ TcBR. Since the spin temperature in the H I region at 
radius r < 0(0.1 Mpc) is very close to the Hi gas tempera- 
ture (Madau et al. 1997), a situation of Ts > Tbg is realized 
if the UV radiation field is not severely shielded. From this 
equation, we obtain 

3 T. 



Ann 



(87) 



leTcBR ' 

On the other hand, if the UV radiation field is shielded, an 
environment of Ts < Tbg leads to 



A 1 3 T. 

" 16 Ts ^ 



(88) 



The line in Fig. 17 was drawn assuming equations (85) 
and using the calculated results of the H I abundance and 
the gas temperature as a function of the time and the ra- 
dius. There is an important difference between the standard, 
soft UV and hard UV cases. In the hard UV case the gas 
temperature has been heated up by the UV background. 
The spin temperature of 21 cm line is then larger than the 
CBR temperature, i.e., Tcbr = 2.725(1 -I- z) K. The 21 cm 
signal is, therefore, an emission, and its amplitude is a fac- 
tor ~ Tcbr /7s smaller than that of absorption in cases of 
Tcbr > Is (Scott & Rees 1990; Madau et al. 1997). 

Figure 18 shows the total flux densities emitted through 
the C II and O I fine structure lines and the H I 21 cm line 
in the calculated domain of r ^ 0.2 Mpc for the case of 
5b = lO**. The high density of the matter causes a slow 
expansion of ionized regions. There are larger volume of C ll 
and O I regions which are affected by the UV pumping than 
in the case of 5h = 10^ . The flux densities are thus enhanced 
in this high density case. 

4.4 Effect of Collisional Excitation 

In dense regions, an effect of collisional excitation on spin 
temperatures of fine structure transitions can not be ne- 
glected. For example, our calculations are performed un- 
der the assumption that the number density of hydrogen is 
riH ~ 2 X 10"^((5b/10^)[(l z)/10]^ cm'^ Collisional ex- 
citations and deexcitations of C l and C ll proceed by in- 
teractions of those species mainly with surrounding H I and 
electron. O l atoms can interact with surrounding H I and 
electron to get excited or deexcited (see Figs. 9 and 11). 
The deexcitation rates taken from (HoUenbach & McKee 
1989) are listed in Table 1. Collisions contribute to spin 
temperatures with their efficiency proportional to collisional 
(de) excitation rates, i.e., 

, 3 
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Figure 18. Same as in Fig 17 for the case of (5b = 10^ 
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(89) 
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These rates are less than spontaneous emission rates of the 
fine structure transitions, and possibly of similar magnitudes 
to UV pumping rates [equations (21) and (80)]. If the gas 
density of the observed region were high, kinetic tempera- 
tures of C I, C II and O I can be imprinted on the spin tem- 
peratures through the collisional excitation (Section 2.1.2). 

Signals originating from collisions need to be small in 
order to detect signals from the pure UV pumping related 
to ionizing sources of the universe. The best site to look for 
is a region which meets both of the following two conditions. 

1. a large region of moderate density. It does not emit 
signals of collisional (de)excitations. The signals from this 
region separate in frequency from those from surrounding 
dense regions where spin temperatures include large contri- 
bution of collisions. 

2. a region where ionizing UV photons are effectively 
shielded and non-ionizing UV photons are not shielded and 
exist abundantly. C I, C ll and O I can exist without being 
ionized, and a UV pumping by redshifted UV photons is 
operative there. 



4.5 Observational Constraint on Physical 
Properties 

When one detect a line signal from a region whose physical 
condition is not known in advance from any different obser- 
vation, it is not clear whether the signal originated from a 
UV dominated region or a collision dominated region. Sig- 
nals of different lines are, therefore, necessary to estimate the 
physical condition. A line signal is determined by a physi- 
cal condition of observed region specified by the density, 
kinetic temperature and environmental UV fiux (absolute 
value and spectrum) . Conversely, a signal of a certain metal 
line could provide a constraint on physical parameters to 
satisfy a condition which possibly realizes the detected sig- 
nal. If signals of more than one different lines are detected 
from the same region, respective constraints on the param- 
eters are obtained since a ratio between excitation rates by 
UV photons and collision are different for respective tran- 
sitions corresponding to different lines. Ideally parameters 
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can then be determined and one can estimate whether the 
signal is from UV photons or from collisions (for individual 
lines) . 

We describe an example of how to constrain physical 
parameters of observed regions for demonstration assum- 
ing that line signals of Cll 91 K, Ol 228K and 99K were 
detected. We note that a more precise estimation of line sig- 
nals is highly desirable. Since the precise propagation of UV 
photons which excite C I, C ll and O I is not addressed in 
this study, it should be studied. For the moment, we assume 
that all UV line photons experience effective scattering of 
C 11 and O I although this assumption would not hold for 
low density region (see Section 4.2.2). 

We suppose that three line signals at z — 8.7 are de- 
tected: Tb - TcBR = (1.2 ± 0.4) X 10"^ mK for Cii 91 K, 
~(1.2±0.4) X 10"^ mK for O I 228K and (9.8±3.3) x 10"^ mK 
for O I 99K. Parameters to be constrained are, for example, 
riHi, ycu, yoi, the UV flux at the Lya frequency, i.e.. Fa, 
and the spectral index, i.e., as. The dependence of signals on 
spectral index (or UV color temperature Tuv) is relatively 
small (see Fig. 2-5). We then neglect it and fix the value 
to be as = 3/2. The carbon and oxygen abundances, i.e., 
yCii and yo^, should be estimated anyhow since the values 
and the ratio between them in high redshifts are not know 
precisely. Detailed theoretical studies on the cosmic chem- 
ical evolution may help this estimation (e.g., Kobayashi et 
al. 2007). We assume that the abundances, j/c„ and yoi, are 
the same as those of present ISM. The physical parameters 
in the present setting is then whi, Tgas and Fa- Effects of 
the UV pumping (Section 4.2.2) and the collision (Section 
2.1.2) are calculated. 

Figure 19 shows a parameter region in the 3D space of 
nn, Tgas and Fa in which magnitudes of observed signals 
are reproduced. The marks of * and x correspond to 
parameter sets which fail to predict observed line signals 
except for only one line signal of C ll 91 K, O I 228 K and 
O I 99 K, respectively. The squares and triangles correspond 
to parameter sets which predicts two right signals (C ll 91 
K and O I 99 K) and (O I 228 K and 99K), respectively. The 
circles correspond to parameter sets which rightly predicts 
all three line signals. 
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Figure 19. An example of estimation of physical parameters, 
i.e., the Hi number density (nu), the gas temperature (Tgas) and 
the UV flux at Lyo frequency (Fa)- The marks of -(-, * and x 
correspond to parameter sets which predict line signals consistent 
with only one observed line signal of C II 91 K, O I 228 K and O I 
99 K, respectively. The squares and triangles correspond to pa- 
rameter sets which predicts two right signals (C II 91 K and O I 99 
K) and (Ol 228 K and 99K), respectively. The circles correspond 
to parameter sets which rightly predicts all three line signals. See 
text for the assumed observational signals and an explanation of 
adopted model. 
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4.6 Redshift Dependence of Signals 

Signals are different in magnitude between sources of differ- 
ent redshifts. As seen in equation (71), the line emission is 



proportional to the UV flux 



[cf. equation (50)] 



and the Hubble expansion rate H{z). When one flxes the 
angle 9 between a source and a position r away from it in 
the direction perpendicular to the line of sight, there is a 
relation between the redshift and the position, i.e.. 



<1 



(90) 



Figure 20 shows the differential antenna temperatures 
in units of mK at the fixed angle oi 9 = 5" from the point 
source of the present setting as a function of the redshift 
under the assumption that the observed region is abundant 
in C I, C 11 and O I in drawing respective lines. 

The shape of curves in this figure is understood as fol- 
lows. If a difference between the spin temperature and the 



Figure 20. Differential antenna temperatures (mK) at the fixed 
angle 6 = 5" in the sky from the direction to the point source as 
a function of the redshift. 



CBR temperature is very small, i.e., ATs 
TcBR, then it is approximately given by 

7b — TcBR = ATsr. 



(91) 



In the case of two level states, for example, if the population 
fraction of the ground state is dominant and the contribution 
of the UV pumping to the excitation is small, then the spin 
temperature is given by 

Am go TcBR 



ATs 



■ exp 



( TcBR ) 



(92) 



no gi r. 

The differential abundance of the excited state is propor- 
tional to the Hubble expansion rate, i.e., Ani cx H{t) oc 
(1 + z)^^"^ [equation (71)]. The abundance of the ground 
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state scales as no oc [1 + z)'^, while the optical depth scales 
as r oc {1 + z)^^^ [equation (19)]. Ultimately there is a rough 
scaling of 



n-TcBR oc {1 + zfexp{T,/ [2.725(1 + z) K]}. 



(93) 



At low redshift the CBR temperature is low, and the dif- 
ferential temperature increases with the decreasing redshift. 
At high redshift the differential temperature increases with 
the increasing redshift mainly through the (1 + z)^ factor 
in equation (93) and the fact that the physical scale corre- 
sponding to a fixed angular scale of object at higher redshift 
is smaller in the ACDM model adopted in the present study. 



4.7 Detectability 

Signals through the C I lines are largest as seen in Fig. 16. 
If the signals of UV photons are not contaminated by other 
signals at around the same redshift, we would detect them^. 
The signal predicted in this calculation is ~ 1 mK at 0.05 
Mpc (11" in the case of z = 8.7) if a neutral C I region exists 
at the distance. This is much larger than the background 
noise measured with WMAP in a frequency range of > 
100 GHz of Tbg < 0.02 mK mainly from vibrational dust 
emission (Gold et al. 2010). Figure 21 shows a simulation 
of a signal of UV photons emitted in the early universe at 
z = 8.7 through the C I 38K line. The contour map of the 
differential antenna temperatures in units of mK is drawn 
on the plane of the angle (arcsec) and the velocity shifts (km 
s~^). For this figure it is assumed that a neutral Cl region 
affected by the UV pumping exists in the region of 0.03 Mpc 
^ r ^ 0.06 Mpc from the point source. 

Concerning the detectability, a signal of ^ 1 mK will be 
seen with the ALMA. For example, a 1 mK signal emitted 
through the C I 38 K line at redshift z = 6.3 will be seen by 
the 12m Array using the Receiver Band of No. 3 for observa- 
tions of extended sources at one sigma level if the following 
parameters are selected: effective bandwidth of 5.3 MHz cor- 
responding to the scale of ~ 0.02 Mpc, beamsize of 3"75, and 
exposure time of 7.6 hr. The ALMA Sensitivity Calculator* 
was used for this estimation. Signals from C ll and O I re- 
gions will be detected by future observations to come after 
the ALMA. 



5 CONCLUSIONS 

The reionization history of the universe at early epoch of 
redshift 2; ~ 10 is not know precisely yet. We study signals 
of the cosmological reionization which would have been left 
on the cosmic background radiation (CBR) through a se- 
ries of excitations of fine structure levels of C i, C 11 and O i 
and following emissions (absorption) of line photons. Since 
the reionization of the universe would have proceeded inho- 
mogeneously in space, regions of neutral or low ionization 
states irradiated by non-ionizing ultraviolet (UV) photons 
naturally exist during the reionization epoch. Non-ionizing 



We do not study an effect of IR radiation of QSO itself on 
a total signal. This contribution should be subtracted by some 
means in an analysis of observational data. 

* http: / / www.eso.org/sci / facilities / alma/obsorving/tools / etc/. 
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Figure 21. Contour map of the differential antenna temperature 
(mK) of the C I 38 K line as a function of the angle (arcsec) and 
the velocity shift (km s~^) when a neutral Cl region affected by 
the UV pumping exists in the region of 0.03 Mpc ^ r 0.06 Mpc. 



UV photons interact with C and O to excite them to unsta- 
ble excited states. This excitation followed by decays of the 
excited states leads to an excitation of fine structure levels. 
Since the UV photons as a source of the reionization produce 
signals through fine structure lines in C I, C 11 or O I regions, 
some information on the reionization may be obtained by 
observations for these lines. 

Essentially, one ionizing photon can produce one line 
photon for a fine structure transition by scattering with C 
and O species. Strong signals are, therefore, emitted at loca- 
tions near strong sources of UV emission with large flux of 
UV photon. We then assume a strong point source such as 
quasi-stellar object (QSO) which starts lighting in the early 
universe, and calculate the evolution of ionized region uti- 
lizing a non-equilibrium chemical reaction network. A rough 
picture of ionization is shown by this calculation that C I, 
C 11 and O l regions irradiated by nonionizing UV line pho- 
tons can exist at locations where ionizing UV photons are 
effectively shielded by dense Hi regions. 

There are two classes of UV photons available for the 
UV pumping of fine structure levels: 

1. photons of frequencies just around the fine struc- 
ture lines at the boundary between ionized region and non- 
ionized region for C I, C 11 and O I, 

2. photons which are emitted as more energetic photons 
at point source and redshifted toward the transition energies 
of fine structures in C I, C 11 and O I regions. 

At ionization boundaries, UV photons which can excite 
the fine structure levels are immediately lost since they are 
used by the UV pumping. It is then predicted that very small 
regions emitting line photons energized by UV sources would 
exist at ionization boundaries. However, the signals of the 
latter class of UV photons is expected to be stronger than 
those of the first. 

Outside the ionization boundaries, redshifted UV pho- 
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tons can leave their signatures possibly over wide region. 
The ratios between UV intensities of respective fine struc- 
ture lines are changed in the relaxation process through 
scattering with the C and O species, and are affected by 
the spin states of C and O. Neutral Hi regions near strong 
UV sources emit signals of UV photons available until the 
relaxation is completed. 

Those signals of UV photons could be contaminated 
by signals of collisional excitations if densities in observed 
regions were large. Regions of intermediate densities might 
be good sites to investigate for signals of the reionization 
since an efi'ectivc UV pumping is possible in denser region 
where an effect of collisional excitation is also larger. 

The dependence of magnitudes of signals on the source 
redshift is shown, and the detoctability of such signals is dis- 
cussed taking the Atacama Large Millimeter/submillimeter 
Array (ALMA) as an example. Although the magnitudes of 
signals depend on physical environments of observed points, 
which are roughly described by many fixed parameters in 
this study, detections of fine structure lino of C I, C ll and 
O I might be possible with the ALMA or future projects. 
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APPENDIX A: NUMBER FRACTIONS OF 
PHOTON IONIZING Hi, He I AND Hell 

When a H II region expands, ionization regions of H ll. He ill 
and Hell form in the order of size, smallest first (Madau 
& Meiksin 1994; Madau et al. 1997). In dense regions of 
(5b = 10^ and 10*, however, the order of ionization fronts 
is found to be Hem, Hii and Hen from the treatment in- 
cluding effects of finite optical depths and recombination 
processes in ionized regions. The cross sections for photoion- 
ization reactions are given by 



crH,M = 6.30xlO-i« (^) 



0-Hei(!^) = 7.42 X 10 

( 



0-HeiiM = 1-575 X 10"'* 



(Al) 



(A2) 



(A3) 



where Ei = 13.60 eV, E2 = 24.59 eV and E3 = 54.42 eV 

arc the threshold energies of Hi, He I and Hell, respec- 
tively (Nakamoto ct al. 2001)'-'. 

Ionizing photons of energies £3 ^ E.^ can react with 
Hi, He I and He 11. The probabilities of photons of energy 
E-f to react with H I and He I are 



gi{E^ 



and 



(1 



+e" 



TlHefHei 



(A4) 



(A5) 



where ti{i>) = nHeO'Hen('^)[^f'(^) ~ r'f^'H* ~ ^^1)] is the 
optical depth between the ionization fronts of H I and He 11, 
and Ati is the interval from the time when the light leaves at 
rf"" to that when it arrives at rf^ , and satisfies the relation 
(neglecting the effect of cosmic expansion), i.e, 

^He„ 



(t- Ati)«ri'(t)-cAti 



(A6) 



Similarly, T2 = [nHaHi(i^) + nHeO-Hen (j^)] [rf (t) - rf^{t - 
At2)] is the optical depth between the ionization fronts of 

He I and H I. At'z satisfies 



rfi(i- Ata) ~rf°i(t) -cAf2 



(A7) 



The reaction rate is proportional to the product of the pho- 
ton flux, the number density of target and the cross sec- 
tion, i.e., oc F{E^)nt8,TgctO'{E.f) . Photons emitted at r = 
are absorbed inside the H l ionization front by a factor of 
e~^^ [right hand side (RHS) of equation (A4)]. Similarly, 
the factor, i.e., e~^^e~'^^ in equation (A5), is for absorption 
inside the He I ionization front. Fraction parts in RHSs of 
the equations show the fractions of photons which are used 
for ionizations of Hi and He I, respectively. The first and 
second terms in the square bracket of equation (A4) corre- 



r < r 



Hej 



spond to the fraction for H r ionization in r, 
and rf"^ ^ r, respectively, while the term in equation (A5) 
is for He I ionization which occurs only in rf^' ^ r. 

The total fraction of photons in this energy range to 
react with H i is 



P3(Hl) 



J^JME,)/E,]gi{E,)dE^ 



fg[L,{E^)/E,]dE^ 



(A8) 



where Li, on v with as = 3/2 is assumed in this paper. 
The fraction to react with He I is 



P3(Hei) = ^ 



(A9) 



IZ[ME^)/E^]dE^ ' 

The fraction to react with Hell is P3(Heii) = 1 — P3(Hi) — 

P3(Hci) accordingly. 

Photons of energies E2 ^ E^ < E3 can react with Hi 
and He l. The probability of photons to react with H I is 

riHO-Hi 



gaiE^) = [1 - exp(-T3)]-|-exp(-T3)- 



,(A10) 



where the first term in the right hand side is the probability 
to react inside the He I ionization front, and the second is 
that to react outside, tj, = nHO-Hj{i^)[rf^'' (t) — rf^{t — Afe)] 
is the optical depth between the ionization fronts of He I and 
H I. Ai3 is determined by 

rf'{t-At3)^rf'''{t)-cAt3. 



(All) 

The total fraction of photons in this energy range to react 
with H I is 



P2(Hl) 



J^^'[L,{E-,)/E-,]g3{E-,)dE-, 
J^^[L4E,)/E,] dE, 



(A12) 



The fraction to react with He: is P2(Hei) = 1 — P2(Hi). 

Finally photons of energies Ei ^ E~^ < E2 can ionize 
only H I, and all are used for the H l ionization. 

The fractions of ionizing photons in three energy ranges 

are /i = 0.589 (Ei ^ E^ < E2), f2 = 0.286 (P2 ^ E-y < 
E3) and /a = 0.125 (£3 ^ E~,) under the assumption of 
as — 3/2. The fractions to ionize H I, He I and He 11 are then 
estimated as follows; /hj = /i + /2P2(Hi) -|-/3P3(Hi), /hcj = 
/2P2(Hei)+/3P3(Hei) and = /3P3(Heii). Generally the 
fraction is given by /, = /fcPfc(i) with k the index for 

energy ranges. 



^ Note that tlieir equation (B5) contains minor typos of wrong 
signs in indexes. 



This paper has been typeset from a T^jX/ IM^^ file prepared 
by the author. 
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